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Abstract. We present results of a Monte Carlo study of the equilibrium dynamics 
of the one dimensional long-range Ising spin glass model. By tuning a parameter ct, 
this model interpolates between the mean field Sherrington-Kirkpatrick model and a 
proxy of the finite dimensional Edward-Anderson model. Activated scaling fits for the 
behavior of the relaxation time r as a function of the number of spins N (Namely 
ln(r) oc N'^) give values of tjj that are not stable against inclusion of subleading 
corrections. Critical scaling (r oc N^) gives more stable fits, at least in the non 
mean held region. We also present results on the scaling of the time decay of the 
critical remanent magnetization of the Sherrington-Kirkpatrick model, a case where 
the simulation can be done with quite large systems and that shows the difficulties in 
obtaining precise values for dynamical exponents in spin glass models. 
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1. Introduction 

The low temperature phase of spin glasses is characterized by an extremely slow 
dynamics. A popular method to study this dynamics consists of Monte Carlo simulations 
of small systems, followed by a finite size scaling analysis, leading to the determination 
of scaling exponents. In this paper we will show the difficulties in obtaining precise 
values for these dynamical exponents with two examples. 

The hrst problem is the equilibrium dynamics of the one-dimensional long-range 
Ising spin glass (IdLR) model. The one-dimensional long-range Ising spin glass model |T] 
with N sites is a generalization of the Sherrington-Kirkpatrick spin-glass (SK) model [2]. 
The spins (cTj = ±1 with equal probability) are placed equidistantly on a circle of 
circumference N, and the Hamiltonian is given by the expression 

n = - ^ , (1) 

where the Jij are independent quenched random couplings distributed according to 

Ji,j = CN{a)eij/rij , = A^/7rsin((i - j)7r/A^) , (2) 

where the Cij are iid random variables with zero mean and unit variance (namely 
E{eij) = 0 and E{eE) = 1). The constant is determined by the normalization 
condition = 1- Depending on the value of the real parameter a > 0, the 

model behaves as follows: 

• 0<cr<l/2. The static properties of the model do not depend on a P 0] and are 
thus the same as for the SK model (obtained for a = 0). In particular the critical 
temperature is Tc = 1 independent of a. This is the SK like region; 

• 1/2 < (T < 2/3. The model is in a mean field phase below some > 0; 

• 2/3 < a < 1. The model is in a non mean held phase below some > 0; 

• 1 < a. The model has no phase transition, namely = 0. 

This model has been extensively studied numerically (both in this version [5]-[8] and 
in dilute versions [9l[in] that are more suitable for efficient numerical simulations) since 
its behavior as a function of a is analogous to the behavior of the Edwards Anderson 
Ising (EAI) spin glass model (the nature of the low temperature of this model has been 
the subject of a controversy for many years) as a function of the dimension d: in the 
IdLR model the value a = 2/3 plays the role of the upper critical dimension of the EAI 
model {ducd = 6 for this model) and the value <7 = 1 plays the role of the lower critical 
dimension {died = 2 for this model). The hope has been expressed that there is an exact 
correspondence between a and d and that there is in particular two values of a dehning 
IdLR models with the same universal properties as the EAI model in dimensions 4 and 
3 respectively. This hope is not substantiated by recent numerical results (using the 
diluted version of the model) however [10] . 
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In a recent paper, Monthus and Garel [a give numerical evidence that the 
equilibrium relaxation time r of the IdLR model with Gaussian distributed quenched 
couplings scales like ii^(ln(r)) oc N'^ with tjj ^ 1/3, like in the SK model (see 
mm and references therein), for all values of cr G [0,1]. This is a surprising result 
since the static properties of the model do depend on a for a > 2/3, and it contradicts 
the later analytical prediction that [16] ip = 1 — a m this region. The results of [ni 
have been obtained for the temperature T = Tc/2, but it is natural to expect that ip is 
temperature independent below Tc. In m the relaxation time is dehned from the long 
time exponential decay of the equilibrium spin spin autocorrelation function qj{t) and 
is obtained numerically, for a given disorder sample J, using an eigenvalue technique 
introduced in [T5|. The systems studied in |TT| are however very small with 6 < N < 20. 
Our purpose is to investigate the question further by direct Monte Garlo simulation of 
the dynamics of the model, a method that allows using quite larger system sizes. It 
is however subject to thermal errors and implies the use of the median to analyze the 
data, instead of the average, due to the presence of a tail of very slow disorder samples, 
whose relaxation time cannot be measured practically with our method. 

The second problem is the critical relaxation of the SK model. This is a particularly 
interesting phenomenon from a numerical point of view since the scaling exponents are 
exactly known and, as there is neither the need to equilibrate the system nor to sample 
thoroughly the phase space, we can simulate quite large systems, by spin glass standards 
at least. 

If a very strong constant magnetic held is applied to a spin glass, the individual 
spins align along this magnetic held. If afterward this magnetic held is switched oh, 
the spin glass relaxes very slowly towards a state of non-zero remanent magnetization, 
with some excess internal energy relative to the internal energy at equilibrium. This 
phenomenon has been the object of detailed experimental studies Ha- 

At the critical point this phenomenon is simpler and well understood 
analytically [HI [T9| for the SK model: Starting at time t = 0 from a conhguration 
where all spins are aligned, both the magnetization m{t), the overlap between two 
clones q{t) and the internal energy e{t) relax algebraically towards their equilibrium 
values with simple exponents, namely m{t) oc q{t) oc t~^‘^, e{t) — e(oo) oc with 
he = dq = 1 and 6m = 5/4. It has been argued that below Tc, m(t), q(t) and e(t) relax 
in two steps, hrst algebraically towards a non equilibrium N dependent value and then, 
on some exponentially large time scale (out of reach of numerical simulations), towards 
their equilibrium value [2011181, however the numerical situation is not clear [2T] . 

The outlook of this paper is as follows: in a hrst section we study the equilibrium 
dynamics of the IdLR model. The precise determination of the scaling law governing 
the growth with N of the relaxation time of this model turns out to be quite difficult. In 
a second section we turn to another problem, namely the decay of the critical remanent 
magnetization of the SK model. This is an enlightening case where the exact values 
of the scaling exponents are known and we can test the precision of the Monte Garlo 
method. In a hnal section we present our conclusions. 
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We measure the value of the spin spin autocorrelation function qj{t) at equilibrium, 
and dehne the relaxation time tj as the (unique) solution of the equation qj{Tj) = 
1/2a/< q^ >j, where the average overlap squared is measured in the same disorder 
conhguration. We have shown [T3j that the alternative dehnition involving the disorder 
averaged < >, namely gj(Tj) = 1/2 a/E(< gives very similar results for the 

system sizes considered. As discussed in [13] the method we are using assumes that 
q{t)/ a/< g^ > ~ Git/r) with some function G(-). 

We now proceed to give some technical details of our simulation; The 
autocorrelation function gj(t) is measured for integer values of the argument inside 
a time window of size Wt,n = 6000 Rt,N) where Rt,n is an integer scale factor adjusted 
in such a way that the window width scales roughly like r when the temperature T and 
the number of spins N vary, specihcally we enforce the relation Wt,n ~ 10 median(r). 
For small values of r however, Rt^n sticks to its lowest possible value Rt,n = 1, and 
the bound that we impose is dehnitively not saturated. Inside this time window gj(t) 
is measured for 60 values of the argument (that are multiple of Rt,n) whose logarithms 
are roughly uniformly distributed. In our computer program the relaxation time tj is 
dehned as the smallest integer tup multiple of Rt,n such that qjitup) < l/2-^< g^ >j. As 
a check of our procedure we have also measured tint, the result of a linear interpolation 
between tdown = tup — Rt,n and tup, truncating the result to the lowest integer, namely 
in C language style notations: 


tint = (int) 


tdown 


Qj(tdown^(tup tdown^ 


Qjit, 


down ) 


Qjitup) 


(3) 


The difference between Rp and Unt gives an idea of the systematic errors induced by 
our selection of values of t for the measurement of q.j{t). Needless to say we must make 
such a selection and cannot compute and store qj{t) for all values of t inside the chosen 
window. We note that in contrast to our approach the method of m gives real valued 
relaxation times. 

We take the quenched couplings to be Gaussian distributed, like in [H] (but in 
the cr = 0 case, where we consider binary distributed couplings, in order to compare 
with [TS] E]), and consider Ndis independent disorder samples, with Ndis = 1024 in 
most cases. We first bring the system to equilibrium, using the parallel tempering 
algorithm, that is currently the best existing algorithm for this purpose. The lowest 
temperature is Tc/2 and we have Nt = 16 values of T separated by a fixed interval 
AT = Tc/lO, namely the highest temperature is T = 2Tc, well inside the paramagnetic 
phase where the dynamics is fast. We perform 2 10^ parallel tempering sweeps (a 
parallel tempering sweep consists of — 1 conditional swaps of configurations plus one 
Metropolis system sweep), the second half being used to measure static quantities like 
< g^ > that we will use later. For the largest systems, we have used 31 values of the 
temperature separated by the fixed interval AT = Tc/20, and perform 4 10^ parallel 
tempering sweeps. We have monitored the time spent by a given Markov chain c 
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at each temperature. With the parallel tempering algorithm as N increases, at hxed 
AT, one notices that these {NtY numbers, as measured inside the simulation, spread 
more and more around their common mean. Very soon some becomes zero and the 
algorithm clearly breaks down. In our simulation, out of many thousand chains, no 
chain has a ratio Uc = maxT({n^})/min'r({n^}) > 64, and only 3 a ration above 32, 
a reasonable result according to our experience with this algorithm, if not ideal. This 
equilibrium procedure is done for two independent copies of the system with the same 
disorder samples, two clones. 

Next we measure (for each clone) the autocorrelation function qj{t) along a long 
chain of length Ngweeps of at least 10^ Metropolis sweeps, and always larger than ten 
window length. Measurements are made every 5Rt^n sweeps. We have checked that a 
ten fold increase in the chain length changes very little to the estimated median and to 
the estimated statistical error on this median, indeed in this simulation the main source 
of statistical error is the fluctuations of the disorder and not the thermal noise Haim. 
The data for qjit) are averaged over the two clones. 

We base our statistical analysis of the dynamics of the model on the median of the 
relaxation time distribution. Considering the median rather than the average leads to 
an immense saving in computer time: we only need the value of r for at least the 50% 
faster disorder samples (samples with the smaller values of r), and can forget about the 
slower samples (that we take to have r = +cxo). Our empirical choice of window length 
turns out to be large enough for this aim. 

To compute the median we hrst sort the Ndis data for r, and dehne the median 
as the average of the values for data number Nfiis/‘^ — 1 and Ndis/‘2 (In our simulation 
the number of disorder samples Ndis is always even). For small system sizes usually the 
values of r for the data number Ndis/‘2 — 1 and Ndis/‘2 are the same (remember that our 
r’s are integer valued), and the median is accordingly not a very good measure. With 
the parameters we have chosen for the time window however this never occurs for larger 
systems (that are what matters for our analysis of the scaling of the relaxation time). 
The statistical errors on median(ln(r)) are obtained from a bootstrap analysis. We have 
done the statistical analysis for both tint and tup prescriptions and found results that 
differ by quite less than the estimated statistical error. This check is important as it 
shows that our selection of values of t for the measurements does not produce signihcant 
systematic errors. The plots presented here have been made using tint- 

In order to compare directly with the results of m we have performed runs at four 
values of a, namely a = 0 (SK model), a = 0.25 (inside the SK like region), a = 0.75 
and a = 0.85 both in the non mean held region. 

Let us start with the IdLR a = 0.75 and T = 0.7 case. We have data for systems 
with N = 8, 16, 24,..., 1024. We ht these data, using gnuplot 4.6, to the usual two- 
parameter form median(ln(r)) = CN^. We have monitored the evolution of the htted 
value of -0 and of the reduced chi squared /ndf, as N^in the lowest value of N included 
in the ht is increased. A stable plateau is reached for when Amin = 128, with a 
reasonable reduced chi squared 0.9. The value of for such a ht (128 < A < 1024) 
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a 

p(j) 

Upgrading 

N 

'max 

ln(r) 

0.75 

Ganss 

Random 

1024 

12.0 

0.75 

Ganss 

Systematic 

512 

11.1 

0.25 

Ganss 

Systematic 

1024 

8.5 

0 

±1 

Random 

768 

8.4 

0 

±1 

Systematic 

768 

7.4 


Table 1. Details of the simulation of the IdLR model (with T = 0.7Tc): cr, coupling 
probability distribution, type of upgrading, largest system size simulated, and the 
median of the relaxation time distribution for N = 512. The systematic upgrading 
leads to a faster dynamics. 


is '0 = 0.20, with an estimated statistical error (±0.003) on the third digit (Bnt the 
systematic errors are larger, as we will see later). As an example of the stability of the 
£t we note that a qnite lower cntoff = 32 wonld give -0 = 0.22 (bnt an nnacceptable 
X^/ndf ~ 6). Fignred] shows a plot of ln(median(ln(r))) as a fnnction of ln{N) together 
with the best two-parameter fit. There are sizable deviations of the data from the fitted 
cnrve for low valnes of N (and a fit to the low portion of the cnrve wonld give a mnch 
larger valne of than the valne we obtain). It is natnral to try a £t inclnding some 
snbleading corrections. The data are however not good enongh for a £t to the snm 
of two power laws, and we have tried a fit to a power law pins a constant, namely 
median(ln(r)) = A + CN^. We find a platean for fj, with a reasonable rednced chi 
sqnared ~ 0.3, for A ^min = 32. We will nse in the following the same valnes of A imin ' 
namely 128 and 32 for fits with two parameters and with three parameters respectively. 
The valne of we obtain with this three-parameter fit is "0 = 0.10 ± 0.01. This is many 
standard deviations away form the previons resnlt. 

We have redone the simnlation nsing systematic npgrading for the dynamics rather 
than random npgrading, still nsing the same thermalized initial conhgnrations (bnt 
with a maximal valne of 512 for N). In the former case the spin confignration at 
time t ± 1 is obtained from the spin confignration at time t by applying the (single 
spin) algorithm to every sites snccessively and in a fixed order In the later case the 
(single spin) algorithm is applied to a spin chosen at random, this procednre being 
repeated exactly N times. It is well known that the former does not satisfy detailed 
balance (bnt do preserve the eqnilibrinm state) whereas the later does preserve detailed 
balance. It is nsnally assnmed that both dynamics lead to the same scaling exponents, 
dne to nniversality, and we assnme that this is trne and compare the resnlts from the 
two dynamics, in order to have an idea of systematic errors, obvionsly not inclnded 
in the gnnplot statistical error estimates. We hrst remark that the difference between 
the random and systematic npdating estimates decreases as Nmin is increased, as it 
shonld. A two-parameter fit of the systematic dynamics data (128 < N < 512) gives 
= 0.22 ± 0.1 and a three-parameter fits (32 < N < 512) gives tp = 0-08 0.01. Both 
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results are in reasonable agreement with their random upgrading counterparts, but leave 
us with some unexplained discrepancy between the two-parameter and three-parameter 
hts. 

The fact that we obtain very small values of with the three-parameter £t indicates 
that may be is indeed exactly zero, and that we have critical scaling, namely r oc 
It turns out that critical scaling hts work nicely. A two-parameter critical £t for random 
update gives p = 2.29 ±0.05 with a fair reduced chi squared (~ 1.4) for Amin = 128, and 
a three-parameter fit gives p = 2.24 ±0.03 for Amin = 32 (see Fig. [2]) with a fair reduced 
chi squared (~ 1.5). Systematic upgrade gives p = 2.15 ± 0.04 and p = 2.13 ± 0.02 
respectively. There is a nice overall agreement between the above four estimates of p, 
with systematic errors of few percents. 

The values obtained for the exponent p are not absurd, it has been argued that 
the non equilibrium correlation length of the 3d 122] and 4d [23] EAI model grows 
like f{t) oc where t is the time since the quench, and z{T) ze z{Tc)Tc/T with 

z{Tc) = 6.86 ±0.16 in 3d and z{Tc) ~ 5.4 in Ad. That f{t) oc inhnite volume 

non equilibrium situation means hand-wavily that at time t the length scales below f{t) 
are equilibrated. Crudely speaking it indicates that at equilibrium on a system of size 
L the relaxation time r fulhlls the relation L oc r^AT)^ namely p = z{T)/d, not far 
from the values of the critical scaling hts of our IdLR data (according to [10] the IdLR 
model with a = 0.79 is a proxy to the 4d EAI model, and z{Tc) ~ 5.4 translates into 
p{T/T, = 0.7) = 5.4/(0.7 d) = 1.9). 

We have made simulations for T = 0.5 with the same value of a but we can 
only simulate systems up to A = 192, where already median(ln(r)) = 15.2, namely 
median(r) 410®, compared to median(ln(r)) = 13.8 for A = 1024 in the T = 0.7 
case. Such a value for A is (jugging from the T = 0.7 case) quite too small to obtain 
sensible results for ip. We also made preliminary runs for a = 0.85 with T = 0.7 and 
0.5 with the same no-go conclusion. 

Simulations for lower values of a are easier. We have data for a = 0.25 and 
T = 0.7, namely in a region where the statics is the same as the one of the SK model. 
We use random upgrading with systems up to A = 1024. Two-parameter hts give 
= 0.25 ± 0.01. Three-parameter hts give = 0.19 ± 0.02 (see Fig. [3]). In this case 
there is no satisfactory critical ht: a two-parameter ht (still with A > 128) has a reduced 
chi squared of 5, and a three-parameter ht (A > 32) a reduced chi squared of 5 again. 

In order to compare with existing results for the SK model we have done a 
simulation of this model with binary distributed couplings,using systematic upgrading 
and A up to 768. A two-parameter ht gives fj = 0.30 ± 0.01 (see Fig. HI), in agreement 
with the existing literature [121 US] E]. A three parameter ht gives however a smaller 
value = 0.17 ± 0.03. Repeating the simulation with random upgrading we obtain 
Ip = 0.26 and ip = 0.17 ±0.02 with two-parameter and three-parameter hts respectively. 
In this case again, no satisfactory critical ht is obtained. 

The conclusion of this section is that our data disagree with the claim that '^ = 1/3 
in the spin glass phase for all 0 < cr < 1. We are left with the following possibilities 
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i) Either the activated three-parameter hts are not trustworthy, and V’ do decrease from 
a value close to 1/3 for the SK case to lower values for larger values of a; ii) Or is quite 
smaller than previously thought, possibly exactly zero in the non mean held region; iii) 
Another possibility is that the coefficient in front of is extremely small as found 
in m for the SK model at T = Tc/20. In order to distinguish between theses scenarios 
one would need some new ideas, as brute force is not a possibility here. 
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Figure 1. Data for ln(inedian(ln(T))) as a function of In(A^) for the IdLR model with 
a = 0.75, T = 0.7 Tc and values of N between 8 and 1024. The quenched random 
couplings are Gaussian distributed and the random updating scheme is used.The green 
line (color on-line) is a two-parameter activated fit of the data for N > 128. 
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Figure 2. Data for median(ln(T)) as a function of In(A^) for the IdLR model with 
a = 0.75 and T = 0.7 Tc as in Fig.[Tl The green line (color on-line) is a three-parameter 
critical fit of the data for > 32. 


3. The relaxation of the critical SK model 

We simulate the SK model with binary distributed couplings, using the Heat Bath 
algorithm with random site updating, starting from a configuration where all spins are 
set to one. The number of sites N ranges from N = 1024 to 2^^ = 131072. We perform 
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Figure 3. Same plot as Fig. [2] but in the SK phase. The critical ansatz fails to 
reproduce the data. Here a = 0.25, T = 0.7 Tc and values of N are between 8 
and 1024. The quenched random couplings are Gaussian distributed and the random 
updating scheme is used. The green line (color on-line) is a three-parameter scaling fit 
of the data for N > 32. 
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Figure 4. Data for ln(median(ln(T))) as a function of In(A^) for the IdLR model with 
a — 0 (SK model), T = 0.7 and values of N between 8 and 768. The quenched 
random couplings are binary distributed and the systematic updating scheme is used. 
The green line (color on-line) is a two-parameter activated fit of the data for N > 128. 
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600 Monte Carlo time steps, and average the results over Ndis = 1310720/iV disorder 
samples. This scaling of N^is is such that the estimated statistical errors are roughly N 
independent. 

In hgure Owe show m(t), the magnetization as a function of t. The data show 
no meaningful hnite size dependence, and are in rough agreement with a m{t) oc t~^^^ 
behavior. They show however some bending that makes the precise determination of 
the exponent 5m ambiguous. Data for e{t) and q{t) are in similar agreement with the 
expected exponents 5e = 5q = 1. 

In hgure Owe emphasize the deviations of the behavior of m{t) from the expected 
law by plotting m{t)t^/^ as a function of t for our largest systems. Note that the ordinate 
range was between 10“^ and 1 in Fig. [5] and is now between 0.6 and 1.5 in Fig. O 
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Figure 5. Decay of the magnetization m{t) as a function of t at the critical point 
for the SK model. The predicted scaling m{t) oc is shown. This is in rough but 

not perfect agreement with the data. 


The ratio clearly depends on t up to t ~ 100, where the noise become overwhelming. 
According to [19], the magnetization should scale with hnite N like 

m{t) = N-^/^ F^it/N^/^) . (4) 

The ratio m{t)t^F hgure [6] shows no hnite size effect but is clearly not independent 
of t. This means that we have sizable scaling violation, and we cannot determine the 
exponent from the small t/N'^F behavior of m{t) in a scaling plot, as is usually done 
(see e.g [25]h 



Figure 6. Ratio of m{t) as in Fig. [5] divided by the expected t behavior, as a 
function of t. 

To go further on, one should study the behavior of m{t) for larger values 
of t This led us to perform additional simulations of the largest systems N = 
16384, 32768, 65536, and 131072 with two improvements: the first is to replace [18] the 
magnetization, namely the overlap between spins at time t and spins at time 0, by the 
overlap between spins at time t and spins at time = 3 (the precise value of t^ being 
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irrelevant), 

m{t) = ^'^ai{t)ai{t^) , (5) 

i 

this should not change the exponents governing the large t behavior, but reduce 
considerably the SN ratio by increasing the signal. Next we average the data for m{t) in 
bins of given integer part of ln(f/4)/(ln(2)/2), this reduce the fluctuations and does not 
affect the large t behavior. We have also multiplied the number of disorder samples by 
16 and the number of time steps by 4 for N = 16384, and 32768. The data can be found 
in Fig. [TJ where both m{t) and t are averages inside bins. The points have been slightly 
shifted (horizontally) in a dependent way, in order to increase the clarity of the hgure. 
Up to t ~ 100 there is no visible hnite size effect (even with the 16 fold increase of the 
number of disorder samples, and the corresponding decrease of the estimated errors) 
and the data are steadily increasing with t. We interpret the data for larger values of 
t as the onset of the asymptotic m{t) oc behavior, and not as a hnite size effect, 
the argument being that if it was a hnite size ehect the data for N = 16384 should be 
below the one at = 32768. The question should hopefully be settled by extending the 
time range to say t = 10000 with excellent precision up to 77 = 131072. This would be 
however a formidable task. Using current data it is interesting to show a scaling plot of 
m{t)t^/^ as a function of t/N'^/^ ( see Fig. [H]). Strong violations of scaling are apparent 
for small values of the argument. The data can be interpreted as a very slow approach 
to scaling as N grows, towards a limit where m{t)t^/^ is independent of t. This shows 
how difficult a precise determination of dm is. Even the modest task of determining 
the exponent governing the hrst non leading term, assuming that the leading term is 
exactly 6m = 5/4 is difficult, as the statistical errors on m{t) increase strongly with t. 
Fig. [8] shows, as an example, the result of a £t of the N = 32768 data to the form 

m{t)t^^* = a — b{N‘^^^/t)^, (6) 

in the range t/N'^^^ > 0.4. Here /i = 0.51 ± 0.01. The estimates of the value of this 
exponent and of the plateau height depend strongly on the range of data included in 
the fit however. 

As a side remark we note that we made some preliminary runs using Metropolis 
instead of Heat Bath and/or systematic updating instead of random updating. It turns 
out that scaling violations (measured for example as the radio of m{t)t^f^ between 
t/N'^/^ = 0.01 and 0.1) are smaller with systematic updating. It is expected that scaling 
violations are not universal and can accordingly depend signihcantly on the choice of a 
dynamics, we have no explanation however why systematic updating gives less scaling 
violations than random updating. 

4. Conclusions 

We study by Monte Carlo method the equilibrium dynamics of the one dimensional 
long-range Ising spin glass model. Varying the parameter a governing the decay with 
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Figure 7. Magnetization m{t) divided by the expected behavior, as a function of 
t for our largest systems. Here m{t) is defined as the overlap between configurations at 
time t and = 3, and the data for both m{t) and t are binned (see text). The points 
have been slightly horizontally shifted in a dependent way, in order to increase the 
clarity of the figure. 



Figure 8. Scaling plot of (as in Fig. [Zl) as a function of The data 

can be interpreted as showing a very slow iV —> oo limiting constant behavior. The 
result of a fit of the N = 32768 data to the form = a — is also 

shown. The data for the largest values of have a negligible weight in the fit. 


distance of the spin spin coupling of this model is similar to varying the dimensionality 
of the canonical Edwards-Anderson Ising spin glass model. We extend up to systems 
with N = 1024 spins the results obtained by Monthus and Garel in [TT], who argued 
that the dynamics is activated with ln(r) oc N'^ with "0 = 1/3 for all values of a, 
from the Sherrington-Kirkpatrick limit {a = 0) to the physical dimension. We hnd 
some unexpected instability in the hts, possibly indicating critical scaling for r (namely 
r oc at least in the non mean-held region. It would mean that the IdLR model 
dynamics is similar to the EAI dynamics in the RSB picture. 

We then turn to an illustrative example of the general difficulties of obtaining precise 
estimates of scaling exponents from numerical simulations of disordered systems, by 
studying the decay of the critical remanent magnetization of the Sherrington-Kirkpatrick 
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model. Here the exact value of the exponent is known and one can simulate systems 
of very large sizes (up to 2^^ here). In this example the asymptotic regime is only 
approached for systems with hundred of thousands spins, that are dehnitively out of 
reach in usual situations, including the equilibrium dynamics of the one dimensional 
long-range Ising spin glass model. 
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